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Abstract 

A systematic approach to Liouville integrable defects is proposed, based on an underlying 
Poisson algebraic structure. The non-linear Schrodinger model in the presence of a single 
particle-like defect is investigated through this algebraic approach. Local integrals of motions 
are constructed as well as the time components of the corresponding Lax pairs. Continuity 
conditions imposed upon the time components of the Lax pair to all orders give rise to sewing 
conditions, which turn out to be compatible with the hierarchy of charges in involution. 
Coincidence of our results with the continuum limit of the discrete expressions obtained in 
earlier works further confirms our approach. 
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1 Introduction 

The issue of integrable defects in discrete and continuum (classical-quantum) integrable sys- 
tems has been the subject of increased research interest during the last two decades or so 
[1]-[TT]. Recently one of us [12] proposed an algebraic approach for the description of a 
Liouville integrable defect in the discrete non-linear Schrodinger model. This approach is 
based on the construction of a A^-site transfer matrix including the defect matrix at a fixed 
point. Classical integrability is guaranteed by the existence of an r- matrix structure 
for the discrete bulk Lax matrices, and the defect matrix. The question of quantum inte- 
grable defects (see e.g. |9l HI] and references therein) will not be treated here, although the 
corresponding formalism is a straightforward variation of the classical one [I2]. Subsequent 
derivations of the Poisson-commuting Hamiltonians and the corresponding time-component 
of the Lax pair following the canonical construction |iT31 E] were given. 

This now leads us to propose a similarly fully algebraic picture for a description of 
a Liouville integrable defect in the continuous non-linear Schrodinger model. We restrict 
ourselves to the case of a single point like defect; extensions of this notion will be commented 
upon in the conclusion section. 

The procedure itself is based on the construction of a suitable continuous transfer matrix 
generating the Poisson-commuting Hamiltonians and their associated time-component V of 
the continuous Lax pair: 
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1. The continuous monodromy matrix is built as a coaction: 



T{A,-A,X)=T+{A,xo,X) L{xo,X) T-{xo,-A,X) 



This of course is the immediate continuum hmit of the discrete defect monodromy matrix 
(see e.g. |12l[T5]). Such monodromy matrices were derived in [16]. The matrices are the 
monodromies of the differential operator d/dx + L{x) where L is the continuous Lax matrix 
L{x) associated to NLS [17], and L is the defect matrix. As in the discrete case, Liouville 
integrability follows from asking that L obeys a quadratic Poisson algebra 

fl.2) 



La{X), Lb{lJ,)j= rab(A-/i), La{X)Lb{lJ,) 

with the same r matrix as the bulk monodromy operators, thereby imposing a strong con- 
straint on the Poisson structure of the dynamical variables parametrizing the defect. 



2. The Poisson-commuting hierarchy of Hamiltonians is then obtained from expansion in 
A^^ of the In of the trace of the monodromy matrix (11. ip . Poisson commutation is for- 
mally guaranteed by the underlying quadratic Poisson structure [H], but must be checked 
against possible divergences due to (5-distributions on a support overlapping the defect point. 



3. The time components of the Lax pair are then computed. They are evaluated sepa- 
rately in the right bulk {xq, A) and the left bulk {—A, xq) and on the defect point -from left 
and right. As in the boundary integrable systems [18] it is required that V-^''(xq ) — )■ V*^^^(xo) 
in order to avoid singular contributions from the zero curvature condition for the Lax pair 
U, V: 



U- V' 



U, V 



0, X Xq. (1.3) 



This translates into sewing conditions {C^^} across the defect relating the right and left val- 
ues of the (j — l)th derivatives of the fields by functions of lower derivatives and the defect 
parameters. Sewing conditions are thus understood as necessary conditions to allow identi- 
fication of the Hamiltonian equations of motion deduced either from H^^^ or from the zero 
curvature condition for the Lax pair U, V*-*-*; in other words they act as "regularizations" in 
the canonical [131 HI] procedure yielding V*-*-* and iJ*^*-* through the classical r-matrix. They 
will be shown in our example to be sufficient conditions. 
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4. Consistency of the procedure then requires to make sure that the sub-manifold of the 
sewing conditions {C^!^} is invariant under the Hamiltonian action, which reads as: 

^n^'\ belongs to the ideal generated by (1.4) 

Once this is checked, we are justified in defining our Hamiltonian dynamical system as a 
Liouville-integrable defect in the continuum. 

As it clearly appears from our construction in Point 1, the monodromy matrix with de- 
fect realizes a Hamiltonian formulation for a Backlund (or rather dressing) transformation 
procedure yielding an "integrable" defect in a Lagrangian approach such as described pre- 
viously in |3] and recently in a very explicit way in [19]. Because we are from the very 
beginning in a Hamiltonian framework involving the r-matrix structure and associated con- 
struction of the Lax pair, we are indeed safe in stating that we are establishing a strong 
Liouville-integrability for our defect theory. In addition this provides the suitable basis for a 
quantization procedure. By contrast integrability discussed in jUdS] can be characterized as 
weaker "Lax" -type integrability, in that they only show the existence of modified conserved 
quantities and their invariance under time evolution triggered by (in the Hamiltonian lan- 
guage) the third Hamiltonian. Higher time evolution cannot be discussed in this Lagrangian 
framework hence Liouville integrability can not be proved. 

Exemplification of the fourfold pattern will be given in Section 2 on the example of NLS 
equation. It must be emphasized that our direct construction can also be formulated as a 
continuous limit of the discrete construction given in [12] . We shall comment on this in a final 
section 3: The continuous limit is here formulated as the replacement of the discrete index 
n by the continuous variable x; the introduction of the normalization scale A in the bulk 
matrices exactly as explained in Section 5 of the previous paper; and the introduction of this 
scale A and an overall factor on the defect Lax matrix (disregarding here the parametrization 
used in [12], which might be relevant in the case of non ultra- local or extended defects). 

Finally, a general proof of the consistency of such sewing conditions as obtained in Step 3 
with the Hamiltonian evolutions triggered by a Lax pair formulation will be given in Section 
4. 

The Hamiltonians obtained by this procedure; the Lax time-operators V(x); and the 
sewing conditions are exactly identified with the direct continuous construction described 
in the previous sections. It appears that the "naive" discrete to continuous limit is here 
consistent, which may be related to the ultra- local nature of the considered theory. 
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2 The continuous NLS model with defect 



The four fold pattern described in the introduction is generic. We shall now exemplify it on 
the simple example of a single point-like defect in a continuum field theory -the non-linear 
Schrodinger model- associated to the Yangian classical r-matrix [20j: r(A) = y, P is the 
permutation operator. 

The starting point in our analysis is the derivation of the corresponding monodromy 
matrix: 

T(A,-A,A) = T+{A,xo,X) L{xo,X)T-{xo,-A,X) 

= Pexpj^^ da;U+(a;)| L(a;o,A) Pexpjy ° da;U-(a;)| (2.1) 

and L satisfy the quadratic relation (11. 2p . We consider the following defect operator 

L(x„) = AJ+ff») (2.2) 

\'j{xo) d{xo) J 

The Lax operator for the NLS model is the familiar (see e.g. [HI [2T]): 

where the bulk fields are canonical i.e. 

V^±(x), ij^iy)] = Six - y), {V^^(x), ^±(t/)} = 0. (2.4) 

Due to the fact that the L satisfies the quadratic algebra (11. 2p the elements a, /3, 7, 6 
realize the following Poison bracket structure: 

|a(xo), /3(a;o)| = /3{xo) 
a{xo), 7(0:0) [ = -l{xo) 



{^(xo), 7M} =2a(xo). (2.5) 

The discussion on the continuum limit of the discrete NLS in the section 3 will further 
justify the present analysis. It will become transparent that the results derived directly from 
the continuum monodromy matrix coincide, as one would naturally expect, with the ones 
obtained as continuum limits of the discrete expressions presented in section 3 . In particular, 



it will be clear (see also [151 122] and expression f l3.19p ) that the continuum analogue of the 
discrete monodromy matrix is given by fl2.ip . 

The continuum "bulk" monodromy matrices satisfy the following differential equation 

?^^^ = V-TH^,y-.X, (2.6) 
and the zero curvature condition is then expressed as: 



V^{x,t) -Y^ {x,t)+ m{x,t),Y^{x,t) =0 x^xo (2.7) 

On the defect point in particular the zero curvature condition is formulated as (this will be 
also transparent when discussing the continuum limit of the discrete theory) 

= V+(xo)L(xo) - L(xo)V-(xo) (2.8) 

at 

and describes explicitly the jump occurring across the defect point. This will be a major 
consistency check of the prescription followed here. The time components at the defect 
point will be explicitly derived below together with the bulk quantities and the "defect" 
quantities V^. 

Now the canonical procedure in extracting the local integrals of motion may be directly 
applied. First we consider the following familiar ansatz for the monodromy matrices: 

T±(x, y; A) = (1 + M/±(a;))e^*("'^)(l + W^{y)y^ (2.9) 

Substituting the ansatz (12. 9 p into (12. 6p . and splitting the resulting equation into a diag- 
onal and an off-diagonal part one obtains 

+ W^V^ - VdW^ + W^V^W^ - = 0, 

ax 
dZ 

^ = U, + U±W^±. (2.10) 

Solution of the latter set of equations provides the explicit expressions of the W^, 
matrices. 

More precisely, let us recall the generating function of the local integrals of motion 

G{\)=\n{trT{\)^ (2.11) 
due to the ansatz (12.91) we can substitute the monodromy matrix accordingly and obtain: 
G{\) = Intr [(1 + iy+(L))e^^(^'"»)(l + W+{xo))-^L{x^){l + iy-(a;o))e^"("«'-^»(l + W-{-L))-A , 

(2.12) 
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but due to the choice of Schwartz boundary conditions at x = ±A we conclude: 

g{X) = \ntr U+(A-o)^i + W+{xo))-^L{xo){l + iy-(xo))e^"("0'-^)) 
Let us first evaluate the first couple of W^^^'s through f l2.10p 

-ij^"{x) + \^^{x)\^^^{x) 

^^"{x) - \^^{x)\^^^{x) 



(2.13) 



(2.14) 



Similarly through fl2.10p the diagonal elements Z*^*-* are given as: 



^+(-1) 

^±(3) ^ 



1 M - 

J dx ■ip^{x)ilj'^{x) 



1 M + a;o 



-A- xo^ 



— f dx ip'^{x)ip'^{x) 
J dx ijj^' {x)'ijj^ (x) 

— J dx ip"^ {x)ijj'^' (x) 
' Jdx U^"ix)tP^ix) - |V'^(a;)|^' 



:..15) 



~ J dx yijj^" {x)ip^ (x) — \'ip^{x)\'' 

The expansion of the generating function Q in powers of j provides the local integrals of 
motion of the model under consideration. Note that due to the fact that for A — > oo the 
leading contribution comes from -keep also in mind that A — )■ oo. It is thus clear that 
the generating function of the local integrals of motion becomes: 

g{X) = Z+(A) + Zfi(A) + ln[(l + W+{xo))-'L{xo)il + l^-(a;o))]ii (2.16) 

and the terms Z^^ provide the left and right bulk charges, whereas the third term of the 
expression above gives the defect contribution. More precisely, the first three integrals of 
motions may be expressed as 

^(1) = 



n 



(2) 



dx if) {x)'ijj (x) + dx ijj^ {x)^p^ (x) + a{xo). 

Jr+ 

dx ip~{x)^~ (x) — dx ip^^x)^^^ (x) 



(2.17) 



- ij^{xo)tp^{xo) + ^+(xo)V' (xo) + 7{xo)ip^{xo) + f3{xo)ip (xo) 



a2(xo) 



(2.18) 



+ (^+(xo)^+(xo))' + 7(xo)V^+'(xo) - /3(a:o)^-'(xo) + ^+'(xo)^-(xo) + 
- V^+(xo)V^~'(xo) - a(xo) (^^{xo)tp^{xo) + (3{xo)^~ixo) + 2^}^ (xq)^' (xq)^ . 

(2.19) 

It is clear that by construction for all the above charges we have: 

|?^(*\ ?/(^')} = 0, (2.20) 

the latter commutation relations have also been explicitly checked for the three charges 
fl2.19p . To show the commutativity of the charges we made use of the exchange relations 
dMl), <UM and also 

{^±(xo), e(xo)} =0, e = |a, /3, 7}. (2.21) 
The fields ip^, at the defect are defined point by analytic continuation: 

ij^ix^) ^ V'^(xo), iJ^ix^) ^ ^P^ixo). (2.22) 

Expressions of the time component V of the Lax pair are known (see e.g. [H]). The 
generic expressions for the bulk left and right theory are given as: 

V+(x, A, /i) = r\X)tra (t+(A X, X)rab{X - fi)T+{x, xo, A)L,(xo, A)T-(xo, -A, A)) 

V"(x, A, /i) = r^{\)tra (TaiA Xo, X)La{xo, \)T^{xo, X, \)rabiX - /i)77(x, -A, A) j 

V+(xo, A, /i) = t-\X)tra {t+{A, Xo, A)r„;,(A - fi)La{xo, A)T-(xo, -A, A)) 

V- (xo, A, /i) = (A)tr, (t+ {A, xo, A)L„(xo, A)r„;,(A - /i)T- (xq, - A, A)) . (2.23) 

In the special case where the r matrix is the Yangian solution the expressions above 
become: 

Y+{x,\,fi) = T+{x,xo)L{xo)T-{xo,-A)T+{A,x) 

A — jJ, 

V-(x,A,/i) = T-{x,-A)T+{A,xo)L{xo)T-{xo,x) 

V+(xo,A,/i) = ^—^L{xo)T-{xo,-A)T+{A,xo) 

A ~ jJ, 

Y-ixo,X,^i) = ^-^T-{xo,~A)T+iA,xo)L{xo). (2.24) 
A — jJ 
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The next step is to expand the latter expressions. Notice that special care is taken by 
construction for the defect point, where separate formulaes naturally emerge. Substituting 
the ansatz for the monodromy matrices we can explicitly derive this expansion. Explicit 
expressions for the first three orders are given below: 

V(/x, x) = °j (2.25) 



xo) 



/i^ — ■?/' {x)ip (x) filp (x) + ijj ' (x) 
fiip~ (x) — ^p^' (x) {x)ip~ (x) 

/i^ — ip~^ {x)ip'^ (x) fi'lp^{x) + ip^' (x) 
fiip^ (x) — ijj^' (x) ijj^ {x)ip^ (x) 



ji'^ - yip^{xo) + /3(xo)jV^ (a^o) IJ'yP^ixo) + /3(xo)j + f(xo) 
/iV'"(xo) - iP~'{xq) {^^{xq) + /3(xo))^"(xo) 

V - '4^^{.Xq){'iIj-{xq) + i{xq)) /iV'+(Xo) + V'+'(xo) 

+7(xo)) +q{xq) V^+(a;o)(^V^^(xo) +7(a^o))^ 



(2.27) 



where we define 



f(xo) = (xo) - a(xo) (^/3(xo) + 2V^+(xo) 
g(xo) = -V^"'(xo) - a(xo) (7(2:0) + 2^"(xo)) . (2.28) 
Due to continuity requirements at the points Xq , Xq (see also a similar argument in [TH] ) 



I.e. 



V+^'^H^) ^ V+(^)(xo), V-('^)(xo ) ^ V-(^)(xo), ^ Xo (2.29) 



we end up with the following sewing conditions associated to the defect point: 



d}' : 




(xo) 


- 


(a;o) - /3(xo) = 0, 










-r 


(xo) - 7(xo) = 








'(xo) 




'{xq) + a(xo)/3(xo) + 2a{xo)i)^ 


"(xo) = 






'(a^o) 




'{xq) + a(xo)7(xo) + 2a{xo)i)~ 


(xo) = 






"(a^o; 


) - V^' 


^"(xo) + 2^-(xo)'^~(xo) 





-2(/3(xo)^+(xo) + /32(a;o))^-(xo) + (/3(xo)7(xo) - 2a\xQ))xlj+ {xo) - [3{xQ)a\xo) 
Cf : V""(xo) - V'+"(xo) + 2V'-(xo)'V'"(xo) - ^+(xo)'V5+(xo) + 2a(xo)^+'(xo) 

-2(^7(xo)V^^(xo) + 7^(xo))V^"(xo) + (/3(xo)7(xo) - 2a^(xo)) V'^(xo) - 7(xo)a^(xo) 

(2 

Higher (j-th) sewing conditions involving jumps of higher (j — l)-th derivatives of the fields 
will arise from the construction of time components of higher Lax pairs. 

Step 4 of the procedure now follows: we need to explicitly check the compatibility of 
the sewing conditions with the hierarchy of Hamiltonian evolutions i.e we show that generic 
relationships of the type (11 ■4p can be implemented consistently with the commuting time 
evolutions. We have in particular checked that: 

{n^'K } = ci^)±c«5'(o), 

+ (^+(xo) + V^-(xo) + 2(3{xo))^^J-{xo)C^y 

(2.31) 

In section 4 we shall formally prove the compatibility of the generic sewing constraints, 
emerging from continuity conditions imposed on the time components of the Lax pairs, with 
the charges in involution. 
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Having defined both the local integrals of motion, and the corresponding Lax pairs we 
now extract the associated equations of motion. Hamiltonian equations on the one hand 

t) = {H^^\ tl)^{x, t)}, t) = {H^^\ ^^(x, t)} x^xo 

e(xo,t) = {H^'\ e(xo,t)}, e G {a, /3, 7} (2-32) 

and the zero curvature conditions (12. 7p . (12.81) on the other hand, give rise as should rightly 
be expected to the same equations of motion. We shall focus on the equations of motion 
emerging from the Hamiltonian "H*^^) (and the Lax pair U^, V^(3)). For the left and right 
bulk theories we obtain the familiar equations of motion from the NLS model 

^±(x,t) 

i^{x,t) 

the dot denotes derivative with respect to time. For the defect point 

d(xo) = 7(xo)V^^'(a;o) + l3{xo)i)^' {xq) - a{xo)'^{xQ)i)'^{xQ) + a{xo)P{xQ)i)~ [xq] 

(5{xo) = 2a^(xo)V'"^(xo) - 2a{xo)ij'^' (xq) + a'^{xo)f3{xo) - f3{xo)j{xo)i''^{xo) 

- f3^{xo)ij~{xo) - 2f3{xo)i'^{xo)i'~{xo) 

7(xo) = -2a{xo)iJ~' (xq) - 2a^(xo)^/'"(xo) - a'^{xo)j{xo) + -f'^ip'^xo) 

+ /3{xoh{xo)^-{xo) + 2-f{xo)^+{xo)^-{xo). (2.34) 

The fact that we end up to the same equations of motion from either the Hamiltonian or 
the Lax pair description for all the points on the line further confirms the consistency of the 
whole process. Indeed, identification of the equations of motion from these two procedures 
confirms that the sewing conditions represent a guarantee that the time-like operators of the 
Lax pair are correctly defined by the Semenov-Tjan-Shanskii expression (I2.23P on the rhs 
and Ihs of the defect point consistently with the hamiltonian evolutions, in other words the 
sewing conditions may be conjectured to represent the necessary and sufficient consistency 
conditions for application of the Semenov-Tjan-Shanskii procedure to inhomogeneous forms 
of monodromy matrices. 



d'^^{x,t) 

dx"^ 
d^ip^ix^t) 

dx'^ 



-2\ij^{x,t)\'i:^{x,t) 
-2\^P^x,t)\'^P^{x,t) 



(2.33) 
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3 Discrete NLS model: the continuum limit 



We shall first briefly recall in this section the main results reported in p!2] • From the discrete 
expressions of the local integrals of motion and the corresponding Lax pairs we shall derive 
a consistent continuum limit which will reproduce the results of the previous section further 
confirming the validity of the proposed process. Let us first recall that the bulk Lax operator 
is given by (see e.g. [T7]): 



LajiX) = \Dj + Aj 

-X, 1 

where = 1 — XjXj and the fields x, X are canonical: 



(3.i; 



jxi, X,}=%. (3.2) 
The defect Lax operator is basically the 5I2 one expressed as: 



■'an 




(3.3) 

the index n simply denotes the position of the defect on the one dimensional spin chain. 
Note that the L matrix is required to obey the same ultra-local Poisson bracket structure 
as the bulk matrices L (13. ip so that integrability is ensured. For this reason the elements 
a, /3, 7, 5 satisfy the following exchange relations: 



ttn, In \ = -In 



7„| = 2a„. (3.4) 

Inserting the defect at the n-th site of the one dimensional lattice the corresponding 
monodromy matrix is expressed as: 

r,(A) = L,;v(A)i:,;v-i(A) . . . LaniX) ■ ■ ■ LaliX). (3.5) 

The trace of the monodromy matrix -the transfer matrix t(A)- as customary provides a 
family of Poisson commuting operators guaranteeing the integrability of the system. The 
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expressions of the discrete integrals of motion obtained in [12] from the expansion of Int(A) 
in powers of j. 



>,2 



1 . q;"^ 
^(2) ^ _ Xj+iXj - 2 5Z ~ Xn+lXn-l - /3„X„_i + 7„x„+i - 

jjtn,n-l jy^n 



3 

(3.6) 



where we define 



-'^n.n-l = -'^n-l — In- (3.7) 

Similarly the time component of the discrete Lax pairs L, were explicitly derived in [12] , 
and the corresponding expressions are recalled below: A^^^ remains the same for all sites, 

(3.8, 



,(2) 



for J 7^ n, n + 1 is given by 



A?V) = ( 1 . (3.9) 



whereas 



U / \7„-A„_i 



Also A^- for j ^ n ± 1, n + 2 is given by 



' -jiXj^i + Xj_iNj_i — Xj_2 "-^j^i-i 
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and 



(3) _ / /i^ + X 



(3) _ I /i^ + Xn+lXn,n-l f-^-Xn+l + a;„+2 " N„,+ iX„_|_i 



—IJ'Xn,n-l — Xn-l + N„_iX„_i + — X„^„_iX„+i 



(3.12) 



where we define 



= - - ttn(7n - 2X„_l). (3.13) 

Notice that the continuum limits of the expressions above 03.131) provide the continuum 
quantities given in (12.281) . 

To obtain the suitable continuum limits of the expressions defined above let us first 
introduce the spacing parameter A in the L-matrix of the discrete NLS model as well as in 
the L matrix of the defect (index- free notation): 

/l + AA-A^xX Ax\ 

L(A) = AA + A|" (3.15) 

Let us now introduce the following notation (see also [121 [151 122]): 

Xj — )■ x~{x), Xj — )■ X~{x), 1 < j < — 1, X E {—A, Sq) 
Xj^x+{x), Xj-^X+{x), n+l<j<N, xe{x+,A), (3.16) 

where xq is the defect position in the continuum theory. Note also that in order to perform 
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the continuum limit we bear in mind that: 

n-l 



^J2f^ ^ r dxf-{x) 



A 5^ /, ^ Idxrix). (3.17) 

j=n+l "^^0 

The continuum hmit of the first integral of motion is then given by f l2.17p . 

Notice that in the first integral we considered terms proportional to A, whereas in the 
second integral the first non trivial contribution to the continuum limit is of order A^. The 
respective continuum quantity reads then as in f l2.18p . The continuum limit of H^^\ after 
taking into account terms of order A^, becomes fl2.19p . It is clear that the expressions 
f l2.17p -( PTT^ were obtained by simply identifying: 

= ip^, = -ij^. (3.18) 

Moreover, in the continuum limit the Lax pair is formulated as: 

Lj-(A) ^I + A U(A, x) + C(A2), Aj^V(x), Aj-+i ^ V(x + A). (3.19) 

The discrete zero curvature condition reads as: 

4(A) = A,+i(A)L,(A) - L,(A)A,(A), j ^ n (3.20) 

which in the continuum limit takes the form (keep terms of order A) [151 122] : 



U- V' + 



U, V 



0. (3.21) 



The Lax pair associated to the first integral is quite trivial f l2.25p . The Lax pairs associated 
to the integrals of motion are derived after taking the following limits: 

Lj L+(x), Af^ Y+^^\x). J e {n + 1, ...N}, xe {x^ , A) 



Lj^L-{x), Af^ ^Y-'^^\x). J G {1, ...n-l}, x e {-A, x, 







L„^L(xo), Ai'^ ^Y-^'\xo), Ci^V+(^-)(xo). (3.22) 

Let us now comment on the zero curvature condition at the defect point. Recall the associ- 
ated discrete zero curvature condition: 



L„(A) = A„+i(A)L„(A) - L„(A)A„(A). 
14 



(3.23) 



The continuum limit of the latter formula, bearing also in mind fl3.22p is given by expression 
fl2.8p . The time component corresponding to "H^^-* (terms of order A) is then given by (12.261) 
And the quantities corresponding to H^^^ (terms of order A^) are given by (I2.27p . The valid 
continuum limits taken above provide extra consistency checks on the results obtained in 
the continuum case in the previous section. 



4 Hamiltonian compatibility for the sewing conditions 

We shall formulate in this section a generic proof on the compatibility of the sewing conditions 
with the time evolutions triggered by the hierarchy of Hamiltonians. 

A formal justification of the closure of sewing conditions on themselves under a linear 
evolution triggered by the integrable Hamiltonians generated by In t{X) can actually be given 
once the linear time evolution of the time-like component of the Lax pair is established. It 
is convenient to start our proof in the frame of discrete integrable models and then consider 
the suitable continuum limit along the lines described in the previous section. 

We recall that the discrete time evolution Aj is defined as follows: given the Lax matrix 
Lj{fi) its time evolution (discrete zero curvature condition) reads as: 

Lj{fi) = Aj+i{X, ^i)Lj{^i) - Lj{ii)Aj{X, /i). (4.1) 

The generating function of the local Hamiltonians (the trace of the monodromy matrix) 
may be expanded as: Int(A) = J2i the generating function A reads as: A_,(A, yu) = 

In the case where the r matrix associated to the system is the Yangian solution the time 
operator may be expressed as (see also [22] and references therein for more details) 

A,(A, ^^) = ^r^r(j - 1, 1; X)T{N,r, X) (4.2) 
X — fi 

where we introduce the notation: 

T(z, j; A) = L,(A)L,_i(A) . . . L,(A), t > j. (4.3) 

Making use of the latter relations f l4.ip -( Hl3l) we deduce the time evolution of A^, (recall 
also {t(A), t(/i)} = 0) i.e.: 



\nt{z), Aj{X,iJ,) \ = ^Aj{z,X), Aj(A,/i) . (4.4) 
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This derivation was given here for simphcity in the discrete framework. It naturally extends, 
-especially given the continuum limit process described in the previous section (see also 
[13 122])^ to the continuous case by: j — )■ x {x = Aj). Equation fl4.4p is valid for all points 
including the defect point n. Similarly its continuum equivalent is valid at every point of 
the interval {—A, A) including the defect point xq. 

Recall now that the sewing conditions are generated by the continuity condition relating 
V^(x^) and V^(a;o). From it becomes possible to write the generic time evolution of 
the sewing conditions: 



Int(z), V±(4,A,/i)-V±(a;o,A,/i) 
V±(x± z,A), V±(xo^,A,/i)] - [v±(xo,z,A), V±(xo,A,/i) 
AV±(z,A), V±(4,A,/x)l + [v±(z,A), AV±(A,/x) 



(4.5) 



where self-explanatorily AY^ = Y^{x^) — Y^{xo). Sewing conditions are obtained from 
expansion of AV^(A, fi) in powers of A~^. The m-th power yields the m-th sewing condition. 
Expansion of the relevant terms in (14.51) yields: 

h 



A;>0 



m>0 



Then equation (14.51) becomes: 



fc>0 m>0 

V^ Cf (A) >r^V±(^-)(4,/i) 



A:>0 



j>0 



\3 



fc>0 



j>0 



A^ 



(4.6) 



(4.7) 



Further expanding quantities C^\X), V^*^-'^(A) in powers of A as 



k-l 



k-l 



V±W(a;o,A) = ^V±('='')(xo)A\ (A) = ^ci''^^ A^ 



i=0 



i=0 



and fixing k and m in (14. 7p leads to the following fundamental relation: 



k-l 



k-l 



i=0 



i=0 



(4.9) 
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Now expanding C^\ij), *)(xq in powers of /i: 

C^\f.) = X^Cf'^V', V^^^n^,/^) = Evi^''Ha;o)/i' (4.10) 

1=0 1=0 

and fixing / in the sums above we conclude 

fe-i fc-i 
I^W^ ^{m,0| ^ ^ y±(m+i,0(a;±)j + J2 [v±('='^)(xo), +''^^] . (4.11) 

i=0 i=0 

C^'''^ are matrices with entries being the constraints of the type fl2.30p or hnear combinations 
thereof. The Poisson bracket of any Hamiltonian "H^*^) with the generic scalar constraint 
is now expressed as finite linear combination of the same scalar constraints, and this 
concludes our formal proof on the Hamiltonian compatibility of the sewing conditions. 



5 Conclusions and perspectives 

Let us now summarize what we have achieved at this time. We have formulated a fully 
Hamiltonian framework for a description of a Liouville-integrable point-like single defect on 
a continuous line in a bulk-integrable field theory. The defect is initially introduced as a set 
of discrete dynamical variables independent of the bulk fields, constrained however by the 
requirement of having a Poisson structure parametrized by the same classical r-matrix as 
the bulk space-like Lax operator. A hierarchy of Poisson-commuting Hamiltonians is then 
derived canonically from the combined bulk-defect monodromy matrix. Sewing conditions 
are then imposed by the requirement that the time-like operators of the Lax pair (describing 
each time evolution associated with each Hamiltonian of the hierarchy) be defined consis- 
tently in the left and right neighborhood of the defect point. They fix relations between 
the defect parameters, and left-right limits at the defect of successive derivatives of the bulk 
fields. 

They were shown, first on the specific example of non-linear Schrodinger equation, then 
following a general algebraic argument based on the r-matrix structure, to be compatible with 
all time evolutions triggered by the hierarchy of Hamiltonians in the sense that the Poisson 
bracket of any hamiltonian with any constraint closes (moreover linearly) on the ideal of 
functions on the phase space generated by the constraints; hence they can be simultaneously 
imposed to the dynamical bulk and defect variables evolving simultaneously under action 
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of the full hierarchy. The hierarchy of Hamiltonians together with the hierarchy of sewing 
constraints thus defines a Liouville-integrable system. 

The sewing conditions, as already emphasized, can be understood as the consistency 
conditions for the existence of a canonical Semenov-Tjan-Shanskii type construction for 
well-defined time-like operators in the Lax pair, i.e. operators such that the zero- curvature 
conditions for the Lax connection give the same equations of motion as the Hamiltonian 
evolution computed from the Poisson structure, both in the bulk and at the defect point. 

Finally as a further consistency check the continuous dynamical equations were identified 
with a suitable scaling limit of the discrete formulation in p^. This first allows to put 
into a clearer perspective a number of previous results on "integrable" defects. As already 
commented upon the Lagrangian approach advocated in many works |l]-[8] has difficulties 
in dealing with the fundamentally Hamiltonian notion of Liouville integrability; it is in a 
sense a "single-time" approach instead of the multi-time approach naturally associated with 
the notion of Hamiltonian hierarchy. "Integrability" in this framework essentially means 
that the constructed conserved quantities are shown to be time-invariant under one single 
time-evolution (the one associated to the first non trivial Hamiltonian of the hierarchy; 
usually the second or third one) which is of course weaker than Liouville integrability. A 
formulation closer to ours can be found in [19j, where the Backlund transformation scheme 
of [8] is rewritten using a combined bulk-defect monodromy matrix similar to ours. A similar 
formulation was proposed earlier by |9j. In both cases however the defect matrix is written 
directly as a function of the limit bulk variables (in our language, this means solving directly 
the sewing conditions to get an on-shell defect matrix). This makes the analysis of the 
Poisson structure (required to speak of integrability) tricky since the on-shell defect matrix 
should now have non-trivial PB's with the left and right bulk monodromy matrices- an 
issue which our "off-shell plus constraints" approach eliminates. Note that by contrast, the 
Hamiltonian formulation of a point like Backlund transformation by Sklyanin [23] is precisely 
of this "off shell" type. 

We can now comment on several possible future developments of our scheme: An ex- 
tension of the point-like single defect approach to other ultra-local integrable field theories 
should not raise too many difficulties at least in its principle. Theories considered in |l]-[8] 
are natural candidates to this extension, which would then clarify the issue of actual Liou- 
ville integrability for the proposed defect models. Multiple point-like defects will a priori 
be described by similar combinations of bulk and defect matrices with independent defect 
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parameters and sewing conditions at each defect point. Extended defects should also be 
considered, in the spirit of the so-called Type II Backlund transformation formalism [6l [19] 
for which one should work at providing a Hamiltonian formulation following the lines of our 
present construction. 

A very challenging question is raised when considering non ultra-local theories. The ultra- 
local form of the Poisson structure in our example considerably simplifies the formulation of 
the time-like Lax operators which are essential to our whole scheme. Non ultra-local PB's 
are naturally [16] associated not with single r-matrices but with r, s pairs parametrizing 
semiclassical reflection algebras. The monodromy matrix structure is more complicated 
(general quadratic form with two matrices) hence the defect Poisson structure will also need 
to be extended; in addition the issue of defect/bulk interaction through crossed Poisson 
brackets (even off-shell!) must be addressed. Construction of the time-like Lax operators 
also becomes then a very non trivial operation (see e.g. |18]). 

Finally, let us comment on possible approaches to quantum integrable defects in the 
continuum. As indicated in the Introduction a construction based on a RTT quantum 
algebra was already proposed in [9], UHl E], and is based on the construction of a quantum 
monodromy matrix on a discrete lattice. Note that the monodromy matrix, and the results 
presented in [12J are apparently valid in the quantum well. 

An alternative approach also exists. It was developed by the Annecy group in general 
cases [23] and particularized to NLS in [25j. It uses an ab initio approach through construc- 
tion of factorizable scattering matrices realizing a "reflection-transmission algebra" . It turns 
out that this approach admits a classical limita and a comparison with our current results 
may be quite illuminating: indeed it is not immediate at this stage how classical sewing con- 
ditions gotten from a Semenov-Tjan-Shanski scheme of Lax pair construction may arise from 
this RT algebra. It is interesting to note that self-adjointness requirements on extensions 
of the quantum NLS Hamiltonian [25] closely resemble the first two continuity conditions of 
the classical time-like Lax operator. 

This opens a new, final avenue of investigation. It is possible to compare this classical 
limit of quantum transmission matrices with classical transmission amplitudes on the classi- 
cal defect. These may be obtained by explicitly solving the bulk-plus-defect NLS equation for 
soliton-like configurations through application of classical direct /inverse scattering methods 



^We are indebted to Eric Ragoucy for pointing out this fact to us. 
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to the Lax pair L and V corresponding to the NLS Hamiltonian with the first two sewing 
conditions. We conjecture that the sewing conditions, being regularity conditions on the 
operator describing the time evolution, will allow for the actual existence of non singular, 
computable classical amplitudes for the soliton moving across to the defect, hence allowing 
for the existence of (at least semi classical) transmission matrices. 
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